Creep motion of an elastic string in a random potential 



O 

o 

(N 

< 
(N 



I 

CZ3 



I 

o 



> 

00 
(N 

oo 

o 

O 



o 



X 
J3 



Alejandro B. Kolton, Alberto Rosso, and Thierry Giamarchi^ 

^ Universite de Geneve, DPMC, 24 Quai Ernest Ansermet, CH-1211 Geneve 4-, Switzerland 

We study the creep motion of an elastic string in a two dimensional pinning landscape by Langevin 
dynamics simulations. We find that the Velocity-Force characteristics are well described by the creep 
formula predicted from phenomenological scaling arguments. We analyze the creep exponent fi, and 
the roughness exponent (. Two regimes are identified: when the temperature is larger than the 
strength of the disorder we find /i ~ 1/4 and ( ~ 2/3, in agreement with the quasi-equilibrium- 
nucleation picture of creep motion; on the contrary, lowering enough the temperature, the values of 
H and ^ increase showing a strong violation of the latter picture. 



Understanding the physical properties of disordered 
elastic systems is a challenging question relevant to a host 
of experimental situations. Indeed such a situation is re- 
alized in many different systems, ranging from periodic 
ones, such as vortex lattices charge density waves 

3, Wigner crystals 0, to interfaces, such as magnetic 
a, U or ferroelectric 0] domain walls, fluid invasion 
in porous media ^3 and domain growth . The com- 
petition between disorder and elasticity in these systems 
leads to unique physical properties and in particular to 
glassy behavior. One particularly important question is 
the response of the system to an external force (magnetic 
or electric field for domain walls, current for vortices, 
etc.). At zero temperature, due to disorder the system is 
pinned and the velocity of the elastic structure remains 
zero up to a critical force Fc- At finite temperature how- 
ever, the barriers to motion due to pinning can always 
be passed by thermal activation an one expects thus a 
finite response at finite force. Although it was initially 
believed that the response was linear [l^l, it was subse- 
quently proposed [l3|, |lj| that due to the glassy nature 
of the disordered system, no linear response would exists. 
The slow dynamics of the system for F <^ so-called 
creep, is expected to be controlled by thermally activated 
jumps of correlated regions over the pinning energy bar- 
riers separating different metastable states. By adding 
some strong assumptions on this physical picture of the 
motion, elegant scaling arguments were used UlEl to 
infer the small F response, leading to the "creep formula" 



V{F) ~ exp 



' T [ F 



(1) 



where Uc is an energy scale, and fi a characteristic expo- 
nent that can be obtained from the relation. 



D-2 + 2C 
2-C 



(2) 



where D is the dimensionality of the elastic system, and 
the exponent ^ is the equilibrium roughness exponent of 
the static system. The above formulas are indeed derived 
under the assumption that the movement is so slow that 
static properties can be used. Relation is remarkable 



since it links the statics with the nonlinear transport of 
a disordered elastic system. 

Going beyond the simple scaling arguments or check- 
ing for such a law has proved to be very challenging. 
Although sub-linear response was clearly seen in various 
systems [ij with good agreement with a precise deter- 
mination of the exponent was clearly more difficult. Rela- 
tion |(2J) has been confirmed experimentally only for mag- 
netic domain walls ^ (see also for vortices). On the 
theoretical side the phenomenological predictions of H1I2|) 
have been derived by a functional renormalization group 
calculation [lil Il9| , starting directly from the equation 
of motion, and valid in an e = 4 — L* expansion. This 
calculation confirmed the phenomenological hypothesis 
made in the scaling derivation and the validity of ^ up 
to the lowest order in e. Although the velocity found 
in the FRG calculation was identical to the one of the 
scaling derivation, important differences were also found, 
notably on the characteristic sizes involved in the motion. 

In spite of these results, the physical picture of creep 
motion is still very phenomenological and many impor- 
tant questions remain open. A systematic study of the 
temperature (or disorder strength) dependence of the 
creep response is particularly lacking, both from experi- 
ments and theory. Moreover, from the theoretical point 
of view, the experimentally very relevant case of low di- 
mensional interfaces (where thermal effects are expected 
to be very important), like elastic lines describing do- 
main walls in thin films, posses a difficult problem to 
tackle analytically, since the FRG [11 113 can hardly be 
used in Z? = 1. Such studies are quite crucial given the 
recent experimental results on creep in magnetic la 13) la 
and ferroelectric [3| systems. Numerical simulations are a 
valuable alternative theoretical tool to address this open 
issue. In this respect, creep simulations of elastic strings 
{D = 1) have been done in the past [21I |2^. but given 
the limited range of velocities available, they were neither 
systematic nor conclusive about the validity of (|1I2|I . 

In this work we used a Langevin dynamics method to 
study the Velocity-Force (V-F) characteristics and the 
dynamic roughness C of an elastic string in a random 
potential. The range of velocities we can explore allows 
a precise check of the creep law. Although we find that 
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the creep law does describe well the data, we also find 
that the equilibrium hypothesis for ^ is not verified. This 
leads at low temperatures or strong disorder to creep 
and roughness exponents that become larger than the 
predicted values = 1/4, ^ = 2/3 respectively. 

We study the creep motion of an elastic string in two 
dimensions driven through a random potential. The 
string is described by a single valued function u{z,t), 
which measures its transverse displacement u from the z 
axis at a given time t. We therefore exclude overhangs 
and pinched-off loops that eventually could be produced 
in the motion of domain walls. Assuming a linear short- 
range elasticity and a purely relaxational dynamics the 
phenomenological Langevin equation describing the mo- 
tion (per unit length) is given by, 

jdtu{z,t) = cd'^u{z,t) + Fp{u, z) + F + ri{z,t) (3) 

where 7 is the friction coefficient, c is the elastic con- 
stant, F is the driving force, and Fp{u, z) = —duU{u,z) 
is the pinning force derived from the disordered pinning 
potential U{u,z). The stochastic force r]{z,t) ensures a 
proper thermal equilibration and satisfies {ri{z,t)) = 0, 
{r]{z, t)ri{z', t')) = 27T(5(z - z')S{t - t'), with (...) denot- 
ing thermal average. The sample to sample fiuctuations 
of the random potential are given by, 

[U{u,z) - U{u',z')]^ = -25{z - z')R{u-u) (4) 

where the over-line denotes average over disorder real- 
izations. In this work we consider a random-bond type 
of disorder, characterized by a short-ranged correlator 
R{u), of range r/ and strength i?(0). A physical realiza- 
tion of this kind of disorder is for instance the random 
anisotropy for magnetic domain walls @ . 

To perform numerical simulations of equation (PJ we 
discretize the string along the z direction, z j = 
0, . . . , L — 1, keeping Uj{t) as a continuous variable. A 
second order stochastic Runge-Kutta method |U 0, 
is used to integrate the resulting equations. To model a 
continuous random potential satisfying Q), we generate, 
for each j, a cubic spline U{uj,i) passing through M 
regularly spaced uncorrelated Gaussian random points, 
with zero mean and variance i?(0)^ 26!]. Moreover the 
random potential satisfies periodic boundary conditions, 

U[u, + M,j) = U{uj,i + L) = U{u„j) (5) 

this defines a finite sample of size {L, M). 

We are interested in the V-F characteristics and the 
dynamic roughness exponent C in the creep regime, for 
different values of T and disorder strength i?(0). The 
average center of mass velocity V is defined as. 



V{F) = 
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FIG. 1: V-F characteristics for R{Q) = 0.30. Curves cor- 
respond to T = 0.24, 0.26, ...,0.42 from bottom (o symbols) 
to top (* symbols). Solid lines are fits of the creep law Q 
with Uc and ^ as fitting parameters. Contrarily to the naive 
creep prediction, the optimal fit parameter ^ is temperature 
dependent. The inset shows log(l/y) vs l3Uc{Fc/FY for all 
T, using their respective fitting parameters /i(T) and UdT). 

The roughness exponent Q is obtained from the average 
structure factor, defined as, 




where q = 2'Kn/L, with n = 1, . . . , L — 1. Prom dimen- 
sional analysis we know that for small q, S{q) ~ q-0-+'^Q 
[27I 1. Fitting our numerical data with this function we 
extract C,. 

We simulate systems of sizes L = 64, 128, . . . , 1024 and 
M = 2L, with c = = 7 = 1. In this paper we show the 
results for L = 512, where finite size effects are negligi- 
ble. We take R{0) = 0.12, 0.30 and temperatures ranging 
from T = 0.8i?(0) to T = 3.5i?(0). We start each sim- 
ulation with a fiat initial configuration (u — 0) at the 
force F ~ Fc/10, and then decrease slowly F in steps of 
AF = O.OlFc up to F/Fc « 0.01. F^ is calculated to 
high precision for each disorder realization using a fast- 
convergent algorithm (26i | . The properties of interest have 
to be calculated when the stationary state is achieved. In 
practice we let the string complete two turns around the 
system for the initial force and one turn for the following 
forces. After this equilibration we estimate numerically 
© and {Tj) approximating the average over disorder and 
thermal realizations by a single time average over one 
turn. 

Typical V-F characteristics obtained in the simulations 
are shown in Fig. ^ In the whole range of temperature 
and pinning strength analyzed we find that the V-F curve 
can be well fitted by the creep formula ^ with Uc and 
fj, as fitting parameters. We thus confirm the predicted 
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FIG. 2: (a) V-F characteristics: (i) R{Q) = 0.12, T = 0.24 (* 
symbols); (ii) R{0) = 0.30, T = 0.30 (o symbols). Solid lines 
are the fitting curves using Q. The inset shows log(l/V) vs 
l3Uc{Fc/FY using the fitting parameters ^ = 0.26, pUc = 3.1 
for (i), and = 0.36, fiUc = 2.2 for (ii). (b) (resp. (c)) 
Structure factor S{q) at F/Fc = 0.02, for (i) (resp. (ii)). 
Sohd lines are fitting curves for small qs. We extract C, ~ 0.67 
for (b) and C, ~ 0.89 for (c). Dashed (dotted) lines correspond 
to the reference value = 2/3 {C,t ~ 1/2). The vertical lines 
indicate the approximate location of the crossover from the 
thermal to the random manifold scaling. 



stretched exponential behavior, being the range of veloc- 
ities in our simulations sufficient to rule out other pro- 
posed forms V2, J8,J. However, contrarily to the naive 
creep relation (0) we find that both fitting parameters, 
and not only Uc as predicted in Refs. HSIlSi can depend 
on temperature. 

Analysis of various values of disorder and temperature 
show essentially two different regimes of creep motion. 
To investigate further these regimes we show in Fig.[21[a) 
the V-F characteristics for two values of disorder and 
temperature representative of each regime. For the small 
disorder case we get the exponent fj, = 0.26±0.01 which is 
compatible with the predicted theoretical value fi — 1/4, 
obtained from ^ using the equilibrium roughening ex- 
ponent Ceq = 2/3 .SO-.Sll . The situation is quite different 
for the strong disorder case, where although the fit with 
the creep formula is still excellent, the value of the 
exponent /i « 0.36 is now clearly in excess with respect 



to the predicted theoretical value. 

To understand in more detail the nature of the two 
regimes we calculate the roughness exponent ( using the 
structure factor O . Quite generally, one can predict that 
the short distance behavior of an elastic string is domi- 
nated by thermal fluctuations (Ct = 1/2)- On the other 
hand, because of the finite velocity, the quenched disorder 
acts effectively as a thermal noise at the largest length 
scale. Thus, in this case, the expected exponent is also 
(^v = 1/2 |33|- Finally, at intermediate length scales, the 
physics is determined by the competition between disor- 
der and elasticity. In particular, in our simulations we 
verified that the Larkin length is negligible. There- 
fore, a random manifold scaling, characterized by a non 
trivial roughness exponent, takes place. In Fig.|2Ib) and 
(c) we show the structure factor for the two cases ana- 
lyzed in Fig. [SJa). As predicted, we get C Ct = 1/2 
for large q. At a certain scale we observe a crossover 
between the thermal and the random manifold scaling. 
The location of this crossover decreases as T (i?(0)) is in- 
creased (decreased) . We can also observe that the second 
velocity-controlled crossover is not achieved in our finite- 
size simulation due to the very slow dynamics. Interest- 
ingly, for the small disorder case, the random manifold 
scaling gives C ~ 0.67, in excellent agreement with the 
equilibrium value Ceq = 2/3, while a much higher rough- 
ness exponent ( « 0.89 is found for the strong disorder 
case. These results are consistent with the previous ones 
for the creep exponent fj,. This conclusion holds for the 
whole range of temperature and disorder strength ana- 
lyzed, as we can see in Fig.l^Ja). We find that the relevant 
parameter to define the two regimes is T/R{0). It would 
be interesting to determine whether such a type of scal- 
ing has a theoretical justification. Moreover, we notice 
that although the values of C and /i depart from the equi- 
librium values, the relation Q seems still to hold, within 
the error bars for the two exponents. This is highly non- 
trivial since equation |5J is derived from a calculation of 
the barriers in an equilibrium situation. 

We discuss finally the temperature dependence of the 
barriers Uc{T), shown in Fig.|S^b). We remark that the 
observed linear temperature dependence is of course pe- 
culiar to the one dimensional wall, where thermal fluc- 
tuations lead to unbounded displacements, contrarily to 
what happens in higher dimensions |20|, |29j . 

In conclusion we have find two regimes of creep motion. 
The first one occurs when the temperature is larger than 
the strength of the disorder, giving ^ ~ 1/4 and ( ^ 2/3 
as predicted by assuming a quasi-equilibrium nucleation 
picture of the creep motion. This implies that the do- 
main wall has time to re-equilibrate between hops, being 
the underlying assumption behind ^ essentially satis- 
fied. The second regime occurs for temperatures smaller 
than the strength of the disorder, and is characterized by 
anomalously large values of both exponents. This clearly 
shows that in this regime the domain wall stays out of 
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FIG. 3: (a) Roughness exponent, C,{!T), and creep exponent, 
/i(T), vs T. C(2^) (* symbols), At(T) (o symbols) correspond to 
i?(0) = 0.30 and C{T) (A symbols), ^i(r) (□ symbols) corre- 
spond to ii(0) = 0.12. The dashed line gives the equilibrium 
roughness exponent C,<tq = 2/3, and the dotted line the purely 
thermal roughness C,t ~ 1/2. The dashed-dotted line corre- 
sponds to the creep exponent predicted from phenomenologi- 
cal scaling arguments, (b) Effective energy barriers Uc(T) vs 
T, for i?(0) = 0.30 (o symbols) and R{Q) = 0.12 (□ symbols). 

equilibrium, and that the naive creep hypothesis does 
not apply. Note that the measured roughness exponent 
is intermediate between the equilibrium value and the 
depinning value C,dep ~ 1-2. The fact that the thermal 
nucleation which is the limiting process in the creep ve- 
locity, is in fact followed by depinning like avalanches was 
noted in the FRG study of the creep . Whether such 
avalanches and the time it would take them to relax to 
equilibrium is at the root of the observed increase of the 
exponent, is clearly an interesting but quite complicated 
open question. 
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